Observation of modified radiative properties of cold atoms 
in vacuum near a dielectric surface 



O 

o 

(N 

< 
o 

(N 



Oh: 



V. V. Ivanov, R. A. Cornelussen, H. B. van Linden van den Heuvell, and R. J. C. Spreeu\\0 

Van der Waals-Zeeman Institute, University of Amsterdani, 
Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 
(Dated: February 2, 2008) 

We have observed a distance-dependent absorption linewidth of cold *^Rb atoms close to a 
dielectric-vacuum interface. This is the first observation of modified radiative properties in vac- 
uum near a dielectric surface. A cloud of cold atoms was created using a magneto-optical trap 
(MOT) and optical molasses cooling. Evanescent waves (EW) were used to observe the behavior 
of the atoms near the surface. We observed an increase of the absorption linewidth with up to 
25% with respect to the free-space value. Approximately half the broadening can be explained by 
cavity-quantum electrodynamics (CQED) as an increase of the natural linewidth and inhomoge- 
neous broadening. The remainder we attribute to local Stark shifts near the surface. By varying 
the characteristic EW length we have observed a distance dependence characteristic for CQED. 

PACS numbers: 42.50.-p, 42.50.Xa, 42.50.Pq 

Keywords: Vacuum field fluctuations, cold atoms, spontaneous emission 
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I. INTRODUCTION 

An electronically excited atom (or molecule) can de- 
cay to the ground state by spontaneous emission. The 
characteristic rate at which this occurs is not simply 
an intrinsic property of the atom but also depends on 
the environment. The spontaneous emission rate is pro- 
portional to the density of electromagnetic field modes 
(DOS, or "density of states"), which is determined by 
the electromagnetic boundary conditions. The DOS can 
thus be modified, and with it the spontaneous emission 
rate. The boundary conditions imposed by the environ- 
ment not only change the radiative linewidth but also 
induce energy level shifts and thus change the transition 
frequencies. These include the electrostatic or Van der 
Waals shift, the Casimir-Polder shift (modification of the 
Lamb shift), and resonant radiative shifts. For a review 
see, e.g. Ref. Q. 

Modified spontaneous emission was first observed by 
Drexhage ^ i2j , using dye monolayers separated from an 
interface by fatty acid layers. Both inhibited and en- 
hanced spontaneous emission have since then been ob- 
served by others in a variety of geometries and circum- 
stances 0,0, EH- Remarkably, the radiative linewidth of 
an atom in vacuo at a distance of the order of an optical 
wavelength from a dielectric surface has never been in- 
vestigated experimentally. Energy level shifts have been 
studied for atom inside cavities 8, 9J and in vapor cells, 
using selective reflection spectroscopy 0, . The situ- 
ation of an atom in front of a distant mirror has recently 
been investigated using a single trapped ion. Both the 
broadening of the radiative linewidth and energy level 
shifts have been reported for this system 0, ■ 
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In this paper we experimentally investigate the radia- 
tive properties of cold (T « 10 //K) atoms of ^^Rb close 
to a glass surface, at a distance on the order of an op- 
tical wavelen gth. Using our method of evanescent-wave 
spectroscopy [T^ "we have observed a broadening of the 
absorption linewidth. We compare our observations with 
calculations based on CQED. 



II. METHOD: EVANESCENT- WAVE 
SPECTROSCOPY 



The radiative linewidth F is proportional to the power 
spectral density of the vacuum field fluctuations at the 
position of the atom i.e. the local DOS. The proxim- 
ity of a dielectric surface imposes a boundary condition 
on the field, changing the DOS. This leads to a modifica- 
tion of F and to energy level shifts 0, [H IH El 111 IH 
IT^ . Both the linewidth broadening and the level shifts 
are significant mainly at distances z < A = A/27r, where 
A is the wavelength of the dominant electronic transi- 
tion. In our case this is the D2 resonance line of Rb, and 
A = 124 nm. 

Therefore we probe the cold atoms near the glass sur- 
face using evanescent- wave (EW) spectroscopy J^] . This 
method is selectively sensitive to atoms very close to the 
surface. An EW appears when our probe beam undergoes 
total internal refiection at the glass surface with index of 
refraction n — 1.51, see Fig. ^ The optical field on the 
vacuum side decays exponentially with the distance z to 
the surface, E{z) cx exp(— z/^). Atoms can absorb light 
from the EW, if their distance to the surface is on the or- 
der of the decay length z ^ X. The decay length can 
be adjusted by changing the angle of incidence 9 accord- 
ing to f (6*) = Xin^ s\T? 6 — 1)~^/^. By adjusting 9 we can 
thus vary the distance scale at which the atoms interact 
with the probe light. By increasing 9 further above the 
critical angle 9c = arcsinn"^, the absorption will occur 
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closer to the surface, where T is more strongly modified. 

For atoms in free space, the absorption profile is given 
by a Lorentzian profile, centered at the (free-space) 
atomic transition frequency cjeg = c/A and with a full 
width at half maximum (FWHM) equal to the natu- 
ral linewidth Foo/Stt — 6.07MHz. Both Ueg and Too 
change in the proximity of the surface. Roughly speaking 
these z-dependent Lorcntzians become convoluted with 
the EW energy density U{z) oc exp(— 2z/^). Therefore, 
in our experiment coegiz) also contributes to the observed 
absorption linewidth through inhomogeneous broaden- 
ing. We expect the width of the resulting absorption 
profile to increase with the angle of incidence 9. In the 
experiment we measured this by tuning an EW probe 
laser across the profile and measuring the absorption. 



III. EXPERIMENT 
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FIG. 1: Scheme of the experiment. A weak evanescent- wave 
probe beam is reflected from the prism surface and collected 
on a photodiode. The reference beam has equal power but by- 
passes the vacuum cell. The difference photocurrent of probe 
and reference beams yields the absorption signal, typically 
a fraction of 10~^ — 10~* of the probe. A second reference 
beam (not shown) was used to monitor variations in the probe 
power. We also collect the fluorescence from the MOT to nor- 
malize for shot to shot variations of the number of atoms. 

The major part of our experimental setup has been 
described previously We produced clouds of cold 

*^Rb atoms using magneto-optical trapping inside a ultra 
high vacuum cell (base pressure p ~ 10~^*'mbar). After 
postcooling in optical molasses we ended up with about 
3 X 10'' atoms, at a temperature of 9 /iK. At this temper- 
ature the Doppler width is 90 kHz (FWHM). The coohng 
lasers were then switched off and the atoms fell down to- 
ward the surface of a glass prism, about 3.6 mm below. 
The center of the cloud reached the prism surface and 
the EW spot after 27 ms. Just before hitting the surface, 
the atoms briefly interact with a weak, p-polarized, EW 
probe beam, see Fig. ^ The intensity of the probe was 
kept well below the saturation intensity to avoid power 
broadening. Using 0.35 ^W and a waist of about 1 mm, 
the maximum saturation parameter was s ~ 0.08. 

The probe beam was derived from a home-built diode 
laser system, locked to the F = 2^F' = (1,3) crossover 



resonance in the D2 line of ®^Rb (780 nm). We used 
an acousto-optic modulator (AOM) to shift the probe 
frequency near resonance with the F ^ 2 ^ F' = 3 
transition and to tune it across the resonance. Before 
sending it into the cell, a fraction of the probe beam 
was split off and sent to a photodiode as a reference. 
After total internal reflection on the prism surface the 
probe beam was focused on a second photodiode. The 
photocurrents of the two photodiodes were subtracted to 
obtain our signal, typically a fraction of 10~^ — 10~* of 
the probe. 

The difference photocurrent was amplified by a 
low-noise current amplifier (Femto, LCA-100K-50M, 
50 MV/A transimpedance) and sent through a low-pass 
filter {RC = 1 ms) to further reduce the noise. All pho- 
todiode signals, including a power monitor and a MOT 
fluorescence monitor were acquired using a digital stor- 
age oscilloscope. The latter two signals were used to nor- 
malize the absorption signals for variations in the probe 
power and shot-to-shot variations in the number of cold 
atoms. 

In Fig. 12 we show a typical EW absorption time trace 
with the probe beam tuned near resonance, for an an- 
gle of incidence di = 9c + 0.52°. The signal has been 
averaged 100 x to reduce the noise. Without filtering, 
the absorption signal has a Gaussian shape due to the 
velocity distribution of the falling atoms. However, in 
order to suppress slow drifts in the difference photocur- 
rent we used AC coupling (i.e. a high-pass filter) on the 
oscilloscope. As a result the Gaussian signal has been 
distorted. Furthermore it is superposed on an exponen- 
tially decaying transient originating from switching off 
the MOT/molasses beams. Although we shielded the 
photodiodes from the molasses light as much as possible, 
some light is still detected. Unfortunately the time be- 
tween switching off the lasers and the arrival of the atoms 
at the surface is fixed by gravity. 

The position of the peak corresponds to the fall time 
of the atoms. The width is given by the ratio of the size 
(~ 3 mm) and velocity {^-^ 0.3 m/s) of the atom cloud as 
it reaches the surface. The height of the peak is ~ 2 mV, 
which corresponds to an absorbed power of ~ 80 pW. 
The time-integrated signal amounts to ~ 3 x 10^ ab- 
sorbed photons, or ~ 2 scattered photons per atom in 
the center of the EW. If we tune the probe laser away 
from resonance, or increase the angle of incidence, the 
signal amplitude decreases and the number of scattered 
photons drops to much less than one per atom. Eventu- 
ally the signal disappears in the noise, which is dominated 
by shot noise. 



IV. DATA PROCESSING 

Despite the signal distortion, we can extract the am- 
plitude and width of the original Gaussian, by fitting the 
filtered time trace to an analytical expression. This ex- 
pression involves error functions due to the known step 
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FIG. 2: A typical time-of-flight signal taken with the 
evanescent-wave probe beam (lOOx average). The peak 
around 27 ms is due to the absorption of evanescent probe 
light by cold atoms arriving at the surface. As a result of 
high pass filtering the signal has been distorted and sits on 
top of an exponentially decaying transient (1/e time 26 ms). 



response function of the filter. We took the fitted height 
of the Gaussian as the measure for the amount of absorp- 
tion. The Gaussian width is essentially constant. Thus, 
the height of the peak is proportional to the absorbed 
EW power, which depends on the EW detuning and the 
angle of incidence. 

For a given angle of incidence we measured time traces 
for different detunings of the EW probe. The fitted Gaus- 
sian height as a function of probe detuning yields an ab- 
sorption profile as shown in Fig. |31 From this we ex- 
tracted a Lorentzian line width by fitting a Voigt profile 
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where A is an amplitude and erfc(.) is the complementary 
error function. This Voigt profile is the convolution of 
a Gaussian with a fixed width A /2tt — 1 MHz and a 
Lorentzian with variable width F. 

The fixed Gaussian linewidth accounts for the finite 
spectral width of the probe laser. Our grating stabilized 
diode laser system has a spectral linewidth coniparable 
to the observed atomic linewidth broadenings [23, l2ll |. 
We determined the laser linewidth in a separate exper- 
iment by observing the beat note between two similar 
but independently locked diode lasers on a photodiode. 
The observed decoherence of the beat signal was well de- 
scribed by a Gaussian with a width of 1 MHz. This is the 
linewidth at short (ms) time scales. For longer timescales 
we rely on the feedback loop of the laser locking electron- 
ics. 

The Lorentzian linewidth F contained in the Voigt pro- 
file is a fit parameter. We performed measurements of the 
absorption profile several times for two different angles of 
incidence (9 — 9c = 0.16° and 0.52°). For each absorption 
profile we find one value for F. The results are shown in 
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FIG. 3: Measured and fitted absorption profiles. The AOM 
frequency is the shift imparted to the frequency of the locked 
probe laser to tune it near resonance. Each data point is 
the fitted Gaussian amplitude of a 100 x averaged time trace 
as in Fig. |21 The solid line is the fitted Voigt profile. For 
comparison, the dashed line shows the free-space Lorentzian 
profile. 



Fig. ^ For larger angles, the absorbed power became 
too small compared to the noise, due to the decreasing 
EW volume. The vertical error bars in Fig.0]are entirely 
determined by the scatter of the datapoints as seen in 
Fig. El 

In the limit of large EW decay length, or 9 ^ 9c, we ex- 
pect F to tend to the free-space value Foq. Unfortunately, 
at angles very close to the critical angle, — 9c < 0.05°, 
the probe beam cannot be treated as a plane wave due to 
the finite diffraction angle. In order to avoid this compli- 
cated situation, we performed an independent check by 
measuring the linewidth in free space. The same probe 
laser was used to measure the absorption by the atomic 
cloud while falling, at a height of 2 mm above the surface. 
A short flash of probe light was used to illuminate the 
atoms. The probe beam containing a "shadow" due to 
absorption by the atoms was recorded on a CCD camera. 
In addition to this picture (a), we took two more: a back- 
ground picture with no probe light (b), and a reference 
picture with probe light but no atoms (c). The optical 
density was then obtained as ln[(a — 6)/(c — b)] and was 
measured for different values of the probe detuning. The 
detuning was varied in a similar way as in the EW probe 
measurements. The linewidth was again determined by 
fitting a Voigt profile. We plotted this free-space value 
in Fig. ^ as the datapoint ioi 9 — 9c = 0. Our mea- 
sured free-space linewidth is in good agreement with the 
known value. This shows that there were no unknown 
systematic broadening effects, such as saturation, stray 
magnetic fields, or residual Doppler broadening. 



V. THEORY FOR r(z) 

We will now compare the measured ^-dependence of 
the linewidth to CQED calculations. When the atom 
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FIG. 4: (Color online) Fitted linewidths for varying angle 
of incidence of the evanescent-wave (EW) probe. The data 
points at 9 — 9c = have been measured in free space, instead 
of with an EW probe. The thin lines are the calculated widths 
based on Eq. JSJ, the thick lines based on an integration of the 
optical Bloch equations (see text). The dashed curves (red) 
show the result when level shifts are not taken into account, 
the solid curves (blue) take into account both broadening and 
level shifts. 



approaches the dielectric surface, both the radiative 
linewidth T and the resonance frequency uJeg change in 
a z-dependent way. Although the latter does not change 
the excited state lifetime, it appears as inhomogeneous 
broadening in the experiment, because the evanescent 
wave performs an integration over z. 

The absorption from the EW probe beam can be cal- 
culated by performing a spatial integration of the photon 
scattering rate over the vacuum half space z > 0. The 
photons that are scattered out of the EW by atoms are 
missing from the reflected probe beam so that the reflec- 
tivity drops below unity. This approach works well if the 
absorption is small (<C 1, i.e. no probe depletion). In 
the limit of low saturation the photon scattering rate is 
r(z)s(z)/2, with the saturation parameter given by 



cU{z) 

s^z) + r^z)/4: 
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surface z, in the notation of Ref. Tt^ (note also 
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Here k — 27r/A, and r'P{u) and r^{u) are the Fresnel re- 
flection coefficients for p and s polarization: 
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In our experiment we probe ®^Rb atoms on the tran- 
sition 55*1/2 (i^ = 2) ^ 5P3/2(F' = 3). An atom in the 
excited magnetic hyperfine state \F',mp) — |3,m) can 
decay to the ground state |2,to — q) with q = 0,±1. 
Choosing the quantization axis perpendicular to the sur- 
face, the q — decay channel is governed by r±, the 
q — ±1 channels by r||. The decay rate for a given sub- 
level 1 3, m) is then given by 



liz) = c,„,or_L(z) + (c,„,_i -I- c„,,,i)r||(z), 



(7) 



where 



is shorthand for the square of a Clebsch- 



Gordan coefficient, Cm,q = (2, m — g, 1, g|3, to) . Note 
that this implies that close to the surface the m states 
have different lifetimes Tt^. The different r„(z) curves 
are shown in Fig. [SJ together with r^(z) and r||(z). 
The curve for |to| ~ 3 is not relevant in our experiment 
because our p-polarized probe does not excite these in- 
states. Our ^''Rb atoms are in a random mixture of all 
five |2,to) states. The probe light is linearly polarized, 
perpendicular to the surface, thus exciting q = transi- 
tions. 



VI. ANALYSIS; COMPARISON WITH THEORY 



where U{z) oc exp(— 2z/^) is the EW energy density and 
Jo = 1.6mW/cm^ is the (free space) saturation inten- 
sity. Note that an increase of r(z) not only increases the 
Lorentzian width but also multiplies into the photon scat- 
tering rate, thus increasing the on-resonance rate. This 
effect tends to favor the detection of atoms near the sur- 
face. 

The modification of the radiative linewidth of an atom 
near a plane dielectric surface has been described the- 
oretically in terms of dipole damping rates rj_ and Fn 
The subscripts (_L, ||) refer to dipoles oriented 
perpendicular and parallel to the surface, respectively. 
The dipole damping rates vary with the distance to the 



Integrating the photon scattering rate Fs/2 over all z 
and averaging over the TO-states, we arrive at the absorp- 
tion of the probe; expressed as a fraction: 



AP 
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6l{z) + TlXz)/^ 



dz, (8) 



where p(z) is the atom density. It is z-dependent due 
to the ground-state level shift that accelerates the atoms 
to the surface. This is well approximated by the Van 
der Waals potential — Ca/z"^, resulting in a depletion of 
the density near the surface according to p(z) — po(l + 
(ziv/z)3)-i/2 with zw = {Cl/mghY'^ « 50 nm. Here 
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FIG. 5: Distance dependence of line widtlis (upper) and level 
shifts (lower) of the relevant hyperfine magnetic sublevels. In 
the upper graph, the dashed curves show the dipole damping 
rates, see Eqs. 31 ■ The curves for the F = \mF\ = 3 
states coincide with that for Fn and do not contribute in the 
experiment. The dotted line is the free-space value. In the 
lower graph the level shifts of the F — \mF\ = 3 states do not 
contribute and are again shown as dashed. 



CI = 5.6 X 10 



"^^ J m'^ is the Van der Waals coeflicient 
and mg h is the kinetic energy with which the atoms fall 
onto the surface. 

The z-dependence of the laser detuning Smiz) in 
Eq. ijHJl accounts for the energy level shifts by unequal 
amounts for the ground and excited states. The shift 
of the ground 5Si/2 state is dominated by the Van der 
Waals shift —C^/z^. The shift of the excited 5P3/2 state 
is more complicated, containing also a resonant compo- 
nent with oscillatory z-dependence. We have used ex- 
pressions for the shifts of both the ground and the excited 
states from Ref . ^16a] , using transition line strengths taken 
from Ref. [S^. We have extended the expressions from 
|16| to account for hyperfine structure. Furthermore we 
have multiplied the results by a factor {n^ — l)/{n'^ + 1), 
because our surface is a dielectric instead of a mirror. 
This is known to be correct in the nonretarded limit [2^ , 
which gives the dominant contribution in the experiment. 

It is evident from Eq. (|SJ) that the absorption profile is 
a convolution of Lorentzians with different widths, am- 
plitudes, and central frequencies. The resulting absorp- 
tion profiles are strictly speaking no longer Lorentzian. 
We have numerically calculated the expected absorp- 



tion profiles using Eq. In practice the deviation 

from a Lorentzian is sufficiently small that we can fit 
a Lorentzian to the calculated profiles. The fitted widths 
are plotted in Fig.0] together with the measured widths. 
In the same Figure we also show the result of the cal- 
culation if we do not take the level shifts into account. 
Clearly, the effect of the level shifts on the observed 
linewidth is comparable with the direct broadening ef- 
fect. 

From the Figure we also see that the observed broad- 
ening of up to about 25% is larger than the calculated 
broadening by about a factor two. This cannot be ex- 
plained by the most obvious sources of spurious broaden- 
ing. These include Dopplcr broadening (< 2%), Zeeman 
broadening due to a spurious magnetic field (< 3%), and 
power broadening (< 0.5%). Furthermore these broad- 
ening mechanisms do not show the observed signature of 
increasing with the angle of incidence. A drift of the laser 
frequency can be excluded by the same argument, plus 
the free-space data point. 

A possible mechanism that would have the correct sig- 
nature is transit time broadening. To investigate this we 
numerically integrated the time-dependent optical Bloch 
equations for an atom moving through the EW field. We 
made the approximation that the atom is a two-level sys- 
tem. First, the known ground state level shift was used to 
solve for the accelerated motion z{t) towards the surface. 
This solution was then used to define a time-dependent 
Rabi-frequency Qi{z{t)) = fiio exp(— 2z(i)/^), and sim- 
ilarly for the detuning d{z(t)), and radiative linewidth 
r{z{t)). Using these time dependent parameters we nu- 
merically integrated the optical Bloch equations to obtain 
the time evolution of the Bloch vector (u{t),v{t),w{t)). 
Note that power broadening is naturally included in 
this method. The number of photons scattered by the 
atom on its way down to the surface was obtained as 
J fli{t)v{t)dt [23. Again an averaging over the magnetic 
m-levels was performed. Finally the probe detuning was 
varied and a Lorentzian fit to the obtained absorption 
profile was performed, as before. The results of the Bloch 
equation approach are also shown in Fig. ^ The two 
calculations yield very similar results. This shows that 
transit time broadening does not explain the discrepancy 
between calculations and measurements. 

As a tentative explanation we invoke the presence of 
local Stark shifts caused by charged or polarized particles 
on the surface. Based on a straightforward model calcula- 
tion we find that a surface charge density of 45e/A^ yields 
a 10% linewidth increase. Remarkably, such a charge 
density corresponds to an average distance between the 
charges of order ^ 100 nm, which is just the distance scale 
to which our experiment is very sensitive. These calcu- 
lations only weakly reproduce the angular dependence 
shown by the data. Recently McGuirk et al. have re- 
ported that Rb adsorbed on a Si or Ti surface generates 
local Stark shifts that were measurable as a change in 
the trapping frequency of their magnetic trap j25l |. The 
authors mention that similar effects on a glass surface 
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like ours are very small. However, their experiment mea- 
sured only changes upon depositing clouds of Rb atoms 
on the surface, whereas our experiment is also sensitive 
to statically present adsorbates. Furthermore, there may 
be other charged or polarized adsorbates on the surface. 
For these reasons local Stark shifts do seem to present 
a likely mechanism to explain our results. Our experi- 
ment is also complementary to Ref. 25] in the sense that 
the latter measures a global effect, whereas our experi- 
ment is sensitive only to local variations of the electric 
fields. Unfortunately we have no detailed information 
about possible adsorbates to make a more quantitative 
analysis. 

VII. CONCLUSION 

In conclusion, we have observed a broadening of the 
absorption linewidth of the D2 resonance line of ^^Rb, 



caused by the surface. Part of the broadening can be ex- 
plained as a combined effect of CQED linewidth broad- 
ening and level shifts due to the proximity of a dielec- 
tric surface. The observed broadening of up to 25% 
was about twice that expected from CQED calculations. 
The likely candidate to explain this discrepancy are local 
Stark shifts due to charged or polarized adsorbates on 
the surface. 
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